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Abstract 

It is shown that the eigenvalue problem for the hermitian Wilson-Dirac operator of for 
a uniform magnetic field in two dimensions can be reduced to one- dimensional problem 
described by a relativistic analog of the Harper equation. An explicit formula for the secular 
equations is given in term of a set of polynomials. The spectrum exhibits a fractal structure 
in the infinite volume limit. An exact result concerning the index theorem for the overlap 
Dirac operator is obtained. 
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The effects of magnetic fields on two-dimensional systems have attracted continual inter- 
ests from condensed matter to elementary particle physics. Magnetic fields not only bundle 
energy spectrum but also rearrange the eigenvalues as they continuously changes, giving rise 
to rich band structure. The most impressive is the energy spectrum of tightly bounded elec- 
tron in a uniform magnetic field. It is well-known that the energy spectrum of this system 
has a fractal structure known as the butterfly diagram |TJ . The key ingredient that underlies 
such variety is the topological nature of magnetic fields. 

Relativistic analogue of the Hamiltonian of tightly bounded electron is the hermitian 
Wilson-Dirac operator (HWDO), which attracts renewed interests in the context of overlap 
Dirac operator (ODO) describing chiral gauge theories on the lattice 0. Magnetic field also 
affects the spectrum of HWDO, giving rise to interesting phenomena such as chiral anomaly 
[[|. We also expect the spectrum has a fractal structure similar to the butterfly diagram. 

In this paper we investigate the spectrum of HWDO for a uniform magnetic field in two 
dimensions. We show that the two dimensional system can be converted to one dimensional 
problem for an arbitrary uniform magnetic field and, as the consequence, the spectrum can 
be characterized by a set of polynomials, which enables us not only to understand the fractal 
structure of the spectrum but also to compute the exact index of the ODO. 

The spectral flow of HWDO for the one-parameter family of link variables was inves- 
tigated in Ref. f|] and the connection to chiral anomaly was elucidated. Recently, one of 
the present author reanalyzed the same system P and gave some exact results on the spec- 
trum for a particular set of uniform magnetic fields, for which the index of the ODO can be 
obtained rigorously. The purpose of this paper is to extend them for an arbitrary uniform 
magnetic field. 

We consider a two-dimensional square lattice of unit lattice spacing and of sides L. The 
HWDO H is defined by 

H^x)=a^(2-m)^x)- £ (l^U^x)^ + ft) + ^f±U;{x - fi^x -//))} , (1) 

where (fi = 1,2,3) are the Pauli matrices and fx stands for the unit vector along the 
/x-th axis. The link variables U^x) and the lattice fermion field ip(x) are subject to periodic 
boundary conditions. In the context of ODO m must be chosen to satisfy < m < 2 in order 
for the correct continuum limit to be achieved. We adopt m = 1 unless otherwise specified. 
We are interested in the spectrum of H for the link variables of the form 

Ui(x) = exp 

where t is an arbitrary parameter and x M stands for the periodic lattice coordinates defined 
by x^ = Xfj_ for < x^ < L and x^ + L = x^. For t = Q being an integer, the magnetic field 
F\2{x) = —iUi(x)U 2 (x + V)Ul(x + 2)[/|(x) becomes a constant and the flux per plaquette 
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normalized by 2ir is given by 

Hence Q is nothing but the topological charge 0. The parameter t connects continuously 
the link variables with various uniform magnetic fields belonging to different topological 
sectors classified by the integer topological charge |6|, |j, [I| . 

The characteristics of the spectrum of HWDO found in || can be summarized as; (1) The 
eigenvalues at integral values of t are separated by several gaps and form clusters. (2) For 
noninteger t the H has in general 2L 2 distinct eigenvalues and rearrangements of eigenvalues 
between the clusters occur in a characteristic way as t increases continuously from an integer 
to the next integer. By carefully inspecting the spectral flows one realizes that for special 
integer values oft — r with r being an arbitrary divisor of L each eigenvalue is exactly r-ply 
degenerate. This can be proven by noting that H can be block diagonalized into r 2sL x 2sL 
matrices each describing a one-dimensional lattice system of degrees of freedom 2sL ||. In 
what follows we will generalize the argument of Ref. ]5| concerning the characterization of 
the spectrum of H for special values t = r or t = L 2 /r to arbitrary integer values of t. In 
doing this it is helpful to note the following general facts; (1) The eigenvalues A of H are 
bounded by |A| < |2 — m\ + 2 [fj]. (2) The spectrum of H is periodic in t with a period 
L 2 . (3) The spectrum is an odd functions in t. Hence, it suffices to analyze the eigenvalue 
spectrum for < t < L 2 /2. 

We first analyze the spectrum for t being an integer multiple of L. Let L be a product 
of two positive integers r and s. Then t can be expressed as t = nL 2 /r = nsL, where n is 
some integer coprime with r. Then we have Ui(x) = 1 and U2{x) = exp[2ninxi/r]. Since 
the link variable is independent of X2 and periodic in x\ with a period r, we can simplify the 
eigenvalue problem for H by the following Fourier transformation 

<p{y\ P, Q) = 4 £ £ e- i9l - ipx Mrl + V, x 3 ) , (4) 

sL 1=0 x 2 =0 

where y ranges between and r — 1. The Fourier momenta p and q are given by 

p = k , q=—J. (fc = 0,l,--., J L-l , j = 0,l,-..,s-l) (5) 

Li S 

H is then block-diagonalized into sL 2r x 2r hermitian matrices h(p, q) given by 

., . (B(p,q) C(p,q) \ ... 

where the first (second) row acts on the upper (lower) component of <p(y;p,q). B(p,q) and 
C(p, q) are defined by 

27my\ 




( B (P, q))v,v = -n s v+i,y + 1 1 - cos ( V + — ] (<>,,,,> - n» y ,y>+l > 

(c(p.ff)w = ki,'+ sin fp+^Vs-ki + i- (7) 
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The Syy, is the Kronecker 5-symbol for < y, y' < r — 1 and satisfies the twisted boundary 
conditions 

= e- iq S yfi , tfg = c ^ f0 . (8) 

We thus find that the original two-dimensional system of degrees of freedom 2L 2 is decom- 
posed into sL one-dimensional systems of degrees of freedom 2r, each described by h(p,q). 
The corresponding eigenvalue problem is just a relativistic analog of the Harper equation 

i- 

The eigenvalues of H are determined by the secular equation det(h(p,q) — A) = 0. It 
takes the form 

det(%,g) - A) = /W(A;p) - t^sin 2 ^sin 2 | = , (9) 

where / r (n) (A;p) = A 2r + • • • is a polynomial of order 2r and is defined by 

/W(A;p) = det(%,0)-A) . (10) 

The q dependent term in (|9|) can be easily found from the the explicit form of h(p,q). For 
n and r being coprime h(p, 0) and h(0, 0) are related by an orthogonal transformation and 
rt n) (\;p) then becomes independent of p. We simply write it as /^(A). Explicit forms of 
f^(\) are given in Ref. for r = 1, • • • , 6. 

Later we need to consider (^) for n and r not necessarily being coprime. In general 
fr i^'iP) depends on p and satisfies the factorization property 

rt n) (X;p)= n(/rHA)-^^sin 2 ^sin 2 ^ , (11) 

where n' = n/s' and r' = r/s' with s' being the greatest common divisor of n and r. The 
easiest way to show this is to consider the case of r = L and t = nL 2 /r = n'L 2 /r' = n's'L. 
Since s = 1, the only allowed value of q is 0. The secular equation (H) then becomes 
fr n \^',p) = 0. On the other hand the same set of eigenvalues must be reproduced by the 
expressions (||P with the substitution r, n — > r', n' and q — > 2irj/ s' for j — 0, • • • , s' — 1. The 
identity (|TTD then follows from these. 

We have shown that the eigenvalues of H for t being any integer multiple of L are 
completely characterized by a set of functions /^(A). We now extend the results to an 
arbitrary integer t. Let n and s be positive coprime integers such that t/L = n/s, then we 
may find a positive integer r satisfying t = nr and L = rs. Denoting X2 by two nonnegative 
integers y and I (0 < y < s, < I < r) as x 2 = si + y, we see that the link variables (|2|) are 
independent of I or, equivalently, periodic in y with a period s. We then define new variables 
V(z;q) by 

*M ^EE^M-^!/) > (12) 
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where q, p and z are defined by 

q= , p= h —=- , z = kL + x 1 . (j = 0, • • • ,r - 1, k = 0, • • • ,s - 1) (13) 

r s L 



Since the shift z — > 2; + sL simply corresponds to k — > + s or, equivalently, to p — > p + 27m, 
we may regard \l/(z;g) as functions of z with a period sL. Here a miracle occurs. In 
the new variables ^f(z; q) defined on the one-dimensional periodic lattice, H is again block 
diagonalized into r hermitian 2sL x 2sL matrices of the form (||) and (0). The concrete 
expressions of B and C are obtained from (|?p by making the following substitutions 

q 

y, r, p, q -> 2, sL, -, . (14) 

s 

The secular equations for the eigenvalues of H at t = nr are then given by f^{\\ q/s) = 0. 
Noting the factorization formula fllTD , we find that all the eigenvalues of H at t = nr are 
determined by the following set of secular equations 

^2 (A) = { ~XT si " 2 sin2 4 , (J = 0,l,---,r-l, i = 0,...y-l) (15) 

where n' = n/s' and r' = r/s' with s' being the greatest common divisor of r and n. 

We see from dl5[) that for each j and / there are 2r's 2 distinct eigenvalues and they must 
lie one by one in the 2r's 2 narrow intervals given by the inequality 

< (-l)^ 2 -V s2 - 4 / r ( l2(A) < 1 . (16) 

Each interval contains exactly rs' eigenvalues. Furthermore, the intervals themselves form 
roughly 2r's 2 /n' clusters. This explains the characteristic features of the spectral flows of H 
found by the numerical investigation ||. We can also count the multiplicity of the eigenvalues 
from (0). In particular the multiplicity at t = rn with n and r being coprime is exactly r. 
In Fig. ^ the spectrum of H for uniform magnetic field is shown. One easily see that the 
feather-shape gaps form a fractal pattern. It looks quite different from the butterfly diagram 
[0]. This is due the special choice of m = 1 and the butterfly-like gaps appear for m ^ 1. 
The spectrum for m — 1/2 is shown in Fig. |2] for comparison. 

In order to understand the appearance of the fractal structure of the spectrum we consider 
the infinite volume limit. The magnetic flux a defined by (§) may be an arbitrary real number 
by considering the limit t, L —>■ 00 with a = t/L 2 fixed. In particular the eigenvalues at 
a = n/r with r and n being coprime positive integers form 2r bands given by the inequality 

< (-l) r - 1 2 r - 4 / r (n) W < 1 > ( 17 ) 

where n specifies how these bands cluster each other. Roughly speaking, n if 2n < r or 
r — n if 2n > r stands for the number of near lying bands. An irrational flux is realized as 
an appropriate limit r,n — >• 00. This implies that the finite number of bands for a rational 
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Figure 1: The spectrum of H for uniform magnetic fields (L = 23, m = 1). 



flux split into an infinite number of tiny bands (maybe a Cantor set). Though the spectrum 
appears smoothly varying with the flux due to the low resolution of the plot, such tremendous 
splittings and focusings of the bands take place continually during a small change of the flux. 
As was argued in Ref. 0, it is possible to find the index of the ODO D defined by 

H , , 

D=1+(T 3 - 



/ H 2 

The index of D is given by 

indexD = Tra 3 - ±Dj = -^Tr-|= , (19) 

where Tr implies the sum over the lattice coordinates as well as the trace over the spin indices. 
Since the trace on the RHS of this expression equal the number of positive eigenvalues of 
H minus the number of negative eigenvalues, we can find indexZ) for constant field strength 
configurations by counting the root asymmetry of the secular equations (|I5fl, where by root 
asymmetry of a polynomial equation we mean the number of positive roots minus the number 
of negative roots. In general the origin A = lies out side of the intervals defined by the 
inequalities ( |16D as can be seen from Fig. |l|. If this is satisfied, it is possible to relate index/} 



with the root asymmetry of fj?J(\) = 0. We thus obtain for t = nr = n'r's' 2 

mdexD = —rs'a^l , (20) 
where stands for the root asymmetry of fj^(X) = 0. 
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Figure 2: The spectrum of H for uniform magnetic fields (L = 23, m = 1/2). 



To understand a^ 1 ' the roots of f^(X) = are plotted in Fig. |. As is easily seen, the 
root with minimum absolute value for each r monotonically increases as r and changes the 
sign from minus to plus at r = 4. We thus find = for r < 3 and = 2 for r > A. In 
the case n > 1 the behaviors of are rather complicated for r < An. However, we know 
cr^ 1 ' = for r = 2n, 3n and = 2n for r = An, 5n, • • • from the factorization relation 
(|ll~D. In fact it is possible to show = 2n for sufficiently large r (> An) rigorously by 
directly evaluating the spectral asymmetry of the hermitian matrix (|6|) for p = q = in the 
large r limit ||. We thus obtain for r's 2 > An 1 or equivalently for t = nr < L 2 /A 

index/} = —rs'n' = —rn = V"F 12 (x) , (21) 

2tt x 

where use has been made of (^|). This is the index theorem for the ODO ( |IBD in the abelian 
gauge background in two dimensions [|, |]. 

We have argued that the HWDO can be block-diagonalized for an arbitary uniform mag- 
netic field and the spectrum is described by a relativistic analog of the Harper equation. 
We have found that the polynomials characterizing the spectrum possess a remarkable fac- 
torization property, from which we can understand the fractal nature of the spectrum and 
count the degeneracy of the eigenvalue. The root asymmetry of these polynomials turned 
out to be related with the index of the ODO. We have established the index theorem for the 
uniform magnetic field |-F 12 (x)| < tt/2. The bound is of course not optimal and depends on 
the choice of the parameter m. 
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Figure 3: Roots of / r (1) (A) = for 1 < r < 70. The two curves indicate the known bound 
|A| > ^jl-{2 + V2) sin^. 
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